Oscillatory convection in binary mixtures: thermodiffusion, 
solutal buoyancy, and advection 

D. Jung, P. Matura, and M. Liicke 

Institut fur Theoretische Physik, 

Universitdt des Saarlandes, 
D-66O4I Saarbriicken, Germany 

(Dated: February 8, 2008) 

Abstract 

The role of thermodiffusive generation of concentration fluctuations via the Soret effect, their 
contribution to the buoyancy forces that drive convection, the advective mixing effect of the latter, 
and the diffusive homogenisation are compared and elucidated for oscillatory convection. Numer- 
ically obtained solutions of the field equations in the form of spatially extended relaxed traveling 
waves, of standing waves, and of the transient growth of standing waves and their transition to 
traveling waves are discussed as well as spatially localized convective states of traveling waves that 
are surrounded by the quiescent fluid. 

PACS numbers: 47.20.-k, 47.20.Ky, 47.54.+r, 47.27.Te 
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I. INTRODUCTION 

Convection in binary miscible fluids like ethanol-water, 3 He — 4 He, or various gas mix- 
tures is an example of systems that show rich and interesting pattern formation behavior 
when driven away from equilibrium. It is paradigmatic for problems related to instabilities, 
bifurcations, selforganisation, complex spatiotemporal behavior, and turbulence. Compared 
to convection in one-component fluids the spatiotemporal properties are far more complex. 
The reason is that concentration variations which are generated via thermodiffusion, i.e., 
the Soret effect by externally imposed and by internal temperature gradients influence the 
buoyancy, i.e., the driving force for convective flow. The latter in turn mixes by advec- 
tively redistributing concentration. This nonlinear advection gets in developed convective 
flow typically much larger than the smoothening by linear diffusion — Peclet numbers mea- 
suring the strength of advective concentration transport relative to diffusion are easily of 
the order thousand. Thus, the concentration balance is strongly nonlinear giving rise to 
strong variations of the concentration field and to boundary layer behavior. In contrast to 
that, momentum and heat balances remain weakly nonlinear close to onset as in pure fluids 
implying only smooth and basically harmonic variations of velocity and temperature fields 
as of the critical modes. Without the thermodiffusive Soret coupling between temperature 
and concentration field initial concentration deviations from the mean diffuse away and in- 
fluence no longer the balances of the other fields [such a system displays close to onset only 
the familiar stationary Benard rolls that are called for historical reasons states of stationary 
overturning convection (SOC)]. 

Hence, the feedback interplay between (z) the Soret generated concentration variations, 
(ii) the resulting modified buoyancy, and (Hi) the strongly nonlinear advective transport 
and mixing causes binary mixture convection to be rather complex with respect to its spa- 
tiotemporal properties and its bifurcation behavior. 

Take for example the case of negative Soret coupling, ip < 0, between temperature and 
concentration fields jjj. Then the above described feedback interplay generates oscillations, 
act the buoyancy difference in regions with different concentrations was identified already 
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as the cause for traveling wave convection. 
Oscillatory convection appears in the form of the transient growth of convection at super- 
critical heating, in spatially extended relaxed nonlinear traveling wave (TW) and standing 
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wave (SW) solutions that branch subcritically out of the conductive state via a common Hopf 
bifurcation, and in spatially localized traveling wave (LTW) states. TW and LTW convection 
has been studied experimentally and theoretically in detail 

The transient oscillatory growth of convection was investigated by numerical simulations 

n n 

|15| . Nonlinear SW solutions were obtained only recently 16]. 

Here we compare these different oscillatory solutions and elucidate common features and 
differences. 

In Sec. |n]we describe the system and our numerical methods for investigating it. Sec. IHII 
deals with extended states of TWs, transient growth of SWs, with SW — > TW transitions, 
and with relaxed SWs. Then in Sec. HYI LTWs are discussed. Sec. El contains a short 
conclusion. 



II. SYSTEM 

We consider a binary fluid layer of thickness d at a mean temperature T and with a mean 
concentration C of the lighter component (e.g. ethanol in water). It is confined between two 
solid horizontal plates perpendicular to the gravitational acceleration g. The layer is heated 
from below with a temperature difference AT to the upper plate. The variation of the fluid 
density p due to temperature and concentration variations is governed by the linear thermal 
and solutal expansion coefficients a = — - p ~ and (3 = ~^§^, respectively. Both are positive 
for ethanol- water. The solutal diffusivity of the binary mixture is D, its thermal diffusivity 
is k, and its viscosity is v. 

The vertical thermal diffusion time d 2 /n is used as the time scale of the system and all 
velocities are scaled by n/d. Temperatures are reduced by the vertical temperature difference 
AT across the layer and concentration deviations from the mean concentration by §AT. 

2 

The scale for the pressure is give n by Then, the balance equations for mass , m omentum, 



heat, and concentration 



171 . 1 1 81 ] read in Oberbeck-Boussinesq approximation 



V-u = (2.1) 

d t u = - V (u : u + p - a V : u) + B ; B = a R (5T + 5C)e z (2.2) 

d t 5T = -V Q = -V ■ [u5T - V5T] (2.3) 

d t 5C = -V J = - V • [u5C - LV {5C - if>6T)) . (2.4) 
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Here, ST and 8C denote deviations of the temperature and concentration fields, respectively, 
from their global mean values T and C; Q and J are the associated currents. B is the 
buoyanc, The Dufoor effect fl Q that provides a coup ltag of concentration gradients 
into the heat current Q and a change of the thermal diffusivity is discarded in ()2.Hj) since it 
is relevant only in few binary gas mixtures [2^ and possibly in liquids near the liquid-vapor 
critical point 23| . 

Besides the Rayleigh number R = ^^-AT measuring the thermal driving of the fluid 
three additional numbers enter into the field equations: the Prandtl number a = v/k, the 
Lewis number L = D/k, and the separation ratio il> = — = —StC(1 — C*)f • Here 
kr = TC(1 — C)St is the t her mo diffusion coefficient |l2( and St the Soret coefficient. They 
measure changes of concentration fluctuations due to temperature gradients in the fluid, ip 
characterizes the sign and the strength of the Soret effect. Negative Soret coupling ip (i.e., 
positive St for mixture like ethanol water with positive a and f3) induces concentration 
gradients of the lighter component that are antiparallel to temperature gradients. In this 
situation, the buoyancy induced by solutal changes in density is opposed to the thermal 
buoyancy. Throughout this paper we consider mixtures with L = 0.01, o = 10, and various 
negative values of ip that are easily accessible with ethanol-water experiments. 

When the total buoyancy exceeds a threshold, convection sets in — typically in the 



form of straight rolls. For sufficiently negative ip the primary instability is oscillatory 24]. 
Ignoring field variations along the roll axes we describe here 2D convection in an x-z plane 
perpendicular to the roll axes with a velocity field 

u(x, z, t) = u(x, z, t) e x + w(x, z, t) e z . (2.5) 

This 2D type of convection is commonly enforced experimentally in convection channels of 
small extension in ^/-direction since the rolls are oriented preferentially perpendicular to the 
channel walls j^. 

To find the time-dependent solutions of the partial differential equations describing con- 
vection we performed numerical simulations with a modification of the SOLA code that is 
based on the MAC method [2^, 13]. This is a finite-difference method of second order in 
space formulated on staggered grids for the different fields. The Poisson equation for the 
pressure field that results from taking the divergence of ()2.2j) was solved iteratively using 
the artificial viscosity method (26]. In addition we used multi-mode Galerkin expansion 
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methods |2j 

The boundary conditions for the fields were as follows: we used realistic no slip condi- 
tions for the top and bottom plates at z — ±1/2 and we assumed perfect heat conducting 
plates. Furthermore, impermeability of the horizontal boundaries for the concentration was 
guaranteed by enforcing the vertical concentration current e z • J to vanish at both plates 

e z ■ J = —Ld z (8C — ip5T) (x, z = ±1/2; t) = . (2.6) 

Laterally periodic extended convection structures were simulated by imposing lateral perod- 
icity in x-direction with periodicity length A = 2n/k. Spatially localized convection states 
were obtained in domains of size up to 160. To measure the strength of the thermal driving 
we use the reduced Rayleigh number r = RJ R° c that is scaled by the critical Rayleigh number 
R° c for onset of pure-fluid convection with the critical wave number k° c . The analytical values 
are R° c = 1707.762 and k® = 3.11632. However, to enable comparison with experimental, 
analytical, or numerical results we scale R by the threshold R° c in question, e.g., of the finite 
differences numerical code j^. 

To characterize convection we monitored among others the evolution of the spatial max- 
imum of the vertical velocity field w m3iX , of the oscillation frequency, and of the reduced 
variance 



M = y/(6&)/(6CZ atd ) (2.7) 

of the concentration field. Note that < M < 1 measures to which extent the binary fluid 
is mixed. The better it is mixed the smaller is the spatial mean (5C 2 ) of the square of 
the concentration deviation from its mean - for spatially localized convection states M is 
defined slightly differently, cf. Sec. IIVI The boundary conditions, however, sustain always 
small Soret-induced concentration gradients against the action of advection and diffusion and 
prevents M from vanishing completely. On the other hand, in the conductive state (denoted 
by the subscript "cond" ) with a Soret-induced vertical concentration stratification given by 
5C con d{z) = there is no advective mixing. Its concentration variance, (SC^ ond ) = ip 2 
is largest yielding M cond = 1. 
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III. SPATIALLY EXTENDED OSCILLATORY CONVECTION 

The complexity of (oscillatory) convection in binary mixtures has been the subject of 
many studies j3] that are too numerous to be discussed here. A recent survey on papers 
adressing the role of the concentration field may be found in the introduction of jl^ . 

In this section we compare characteristic properties of spatially extended oscillatory con 



vection in the form of relaxed TWs 



I9I 30, 31 1, SWs |3] and of oscillatory transients [l^ 



into, say, a nonlinear TW. In Fig. ^ we show the bifurcation diagrams of maximal vertical 
flow velocity w max and oscillation frequency uo versus reduced Rayleigh number r for the 
representative Soret coupling if) = —0.25. 

At r osc = 1.335 a SW solution as well as a pair of symmetry degenerate left and right 
traveling wave solutions bifurcate out of the conductive state with wave number k = n. At 
this Rayleigh number the system shows a subcritical Hopf bifurcation with a Hopf frequency 
^H{k = Tr,r osc ) = 11.23. Here two facts are worth mentioning: (i) the bifurcation threshold 
and the frequency are practically the same as for the critical waves with the critical wave 
number k c = 3.135 and (ii) TWs with wavelength A ~ 2 are often observed in large-scale 
experimental setups. This value corresponds to the minimum of r osc (k). Fig. 121 shows 
r osc(k) for different if). The respective onset values for the phase velocities Vp(k) = uj H {k)/k 
depend only weakly on the wave number k. The curves v^{k) are shifted to higher values 
with increasing negative if) due to the higher Soret-induced concentration contrast in the 
conductive state. The saddle node lines rj w (k) are more or less downwards shifted versions 
of r osc (k). The size of the interval, r osc (k) — rj w (k), of subcritical TWs increases when 
if) becomes more negative. For the parameters of Fig. El the TW branches end at large 
amplitudes in a well mixed state of vanishing frequency by merging with the stationary 
SOC solution branch of the same wavelength. For the parameters of Fig. ^ this happens at 
r* = 1.65 where the SOC solution branch is included for convenience as well. 

In contrast to the TW solutions, the SWs remain unstable at least against a transition 
into a (transient) TW. Their solution branches (cf. Sec. lIIID|) end in period dubling cascades 



[16]. 
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A. Advective mixing and diffusion of concentration fluctuations 

When moving along the TW and SW solution branches in Fig. ^ the frequency u and the 
TW phase velocity v p = uo/k decrease monotonically starting from the large Hopf value at 
the bifurcation threshold. This is related to the advective reduction of the Soret generated 
concentration gradients. In fact there is a universal linear relation between the degree of 
mixing and the frequency u which holds for weakly nonlinear as well as for strongly nonlinear 
states and for transients alike. See for example (30 , Q for TWs, Q for LTWs, and Q for 
transients. This scaling relation u / ujh — M holds also for SWs as one can infer from Fig. 01 
Also TWs with wavelength different from 2 show roughly this scaling. 

In particular on the stable, strongly nonlinear TW solution branch the binary fluid gets 
with increasing w max and thermal driving r more and more mixed and M reduces almost 
to zero - only concentration gradients in the narrow boundary layers survive. This holds 
even more so for the well mixed SOC states (see e.g. the SOC2 in Fig. 0] which is marked 
in Fig. by the filled triangle) that resemble closely the corresponding stationary state in 
the pure fluid with the same a and r. 

The transition to convection at r osc is hysteretic, i.e., of first order because the Soret 
coupling coefficient if) = —0.25 is sufficiently negative. The associated precipitous growth of 



19 



convection is caused by an interplay between the solutal contribution to the buoyancy 
that tends to stabilize the conductive state and the effect of advective mixing. The latter 
enhances convection by reducing the adverse effect of the Soret generated concentration 
variations. 

There occurs an "S"-shaped deformation of the unstable part of the TW bifurcation 
branch [dashed (blue) line in Fig. Q when the advection velocity w max has grow n to become 
equal to the TW phase velocity v p . Then the first closed streamlines appear [31 1 in the frame 
of reference that is comoving with the TW phase velocity and where the TW solution is time 
independent. This is indeed a characteristic feature of strongly nonlinear TW convection, 
cf. further below. 

For w max < \v p \, i.e., closer to the threshold r osc all streamlines are open and the field 
structures look almost harmonic - see the TW1 plot in Fig. HJ This weakly nonlinear TW, 
marked by an open circle in Fig. Q is just at the border w max ~ \v p \. This condition marks 
also the border line beyond which a straightforward small-amplitude expansion around the 



February 8, 2008 8 

convective onset breaks down In addition the generic transient growth dynamics of 

oscillatory convection (consisting initially of oppositely traveling waves of roughly equal 
velocity amplitudes \Ar\ ~ |Al|) undergoes a dramatic change that cannot be described at 
a» with ampntude equation mode.s when the now amputudes \A R , L \ approach the border 
line value of u/k ■ 

The qualitative change in the flow topology between weakly and strongly nonlinear TWs 
causes a different mixing behavior with increasing amplitude. This is the reason for the "S"- 
shaped deformation of the unstable part of the TW bifurcation branch [dashed (blue) line 
in Fig. p. It becomes more pronounced with if) becoming more negative (see e.g. Fig. EJ). 
Eventually the strongly nonlinear TW solution branch develops a bistable part there on 
which fast TWs are located 

For relaxed TWs with < \v p \ < w max the remaining open streamlines are spatially 
correlated with inner boundary layers of the concentration field. As indicated in the TW2 of 
Fig. E] they meander along the (green) regions of mean concentration, SC = 0, between and 
around the roll-like regions of closed streamlines. In a right-propagating TW the regions 
of closed streamlines for the right (left) turning fluid domains are rich (poor) in the lighter 
component - here ethanol - and they are displaced towards the upper cold (lower warm) 
plate, where the Soret effect maintains a boundary layer with alcohol surplus (deficiency). 
In addition the meandering open streamlines lie between the upper (lower) closed streamline 
regions and the opposite bottom (top) concentration boundary layer. This structure of 
closed and open streamlines causes the top (bottom) boundary layer to feed high (low) 
concentration only into the right (left) turning roll domain at the location of downflow 
(upflow). Then the fluid becomes diffusively homogenized in the closed streamline regions 
of the rolls leading to anharmonic concentration profiles of trapezoidal shape shown in the 
top plot of the left column of Fig. HJ The motion of the rolls with their specific concentration 
distribution implies a mean concentration current which is directed to the right in the upper 

nnn 

and to the left in the lower half of the layer. See refs. [13, |4£J, yjj for a more detailed 
discussion of this current. 

With increasing w max and decreasing \v p \ the regions of closed streamlines grow at the 
expense of the open ones. Thereby, the former also come closer to the respective opposing 
boundary layer. This decreases the asymmetry of the boundary layer feeding into oppositely 
turning rolls. As a consequence the concentration contrast between adjacent TW rolls 
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decreases until in the SOC state with v p = the rolls are fed symmetrically by both boundary 
layers and mirror symmetry between the rolls is established. 

B. Symmetries 

SOCs and SWs are laterally mirror symmetric around positions of maximal up- and 
downflow, e.g. x=0 in Fig. HJ and the node locations of the fields are fixed in time. This 
symmetry is broken in TWs - see, e.g., the concentration contrast between left and right 
turning rolls in Fig. HI But all fields of SOCs, SWs, and TWs have at every instant definite 
parity under the mirror-glide (MG) operation (x,z) — > (x + A/2, —z) of vertical reflection 
at mid-height, z=0, combined with lateral translation by half a wavelength. We did not 
observe SWs without this symmetry - perturbations breaking it that we introduced for test 
purposes always decayed rapidly to zero. 

Furthermore, all transients investigated in obeyed the MG symmetry with the ex- 
ception of the very early stage in cases where the imposed initial conditions were not MG 
symmetric. But even then the MG symmetry was rapidly restored by a fast decay of MG 
symmetry-breaking modes. Also the transient growth seen in the experiments in an annular 
geometry by Winkler and Kolodner j2| was locally MG symmetric. A time dependent gen- 
eralization of this symmetry was found to be realized in LTW states 143] and an extension to 

n 

3D patterns was observed in |44|. Furthermore, MG symmetric convective structures were 
not only observed in Soret driven convection with only temperature gradients imposed but 
also in thermosolutal convection 

Thus, the MG symmetry that is displayed by the basic conductive state and by the 
linear critical convective modes seems to be quite robust and also persistent in non-linear 
convective structures of pure fluids and of mixtures. The robustness of this symmetry is 
remarkable given that the nonlinearity in the concentration balance, i.e., the Peclet number 
w/L is quite large - of the order of 1000. However, the concentration field is "tamed" by 
being coupled to the velocity and temperature fields. Their shape remains mostly harmonic 
for small supercritical thermal driving like in pure fluid convection. And thus the increase 
in structural complexity associated with a MG symmetry breaking does not occur. 

Finally, the SW fields start at onset to have in addition a definite mirror-timeshift sym- 
metry (MTS), e.g., f(x,z,t) = —f(x,—z,t + r/2) for / = 5C,ST, and the vertical veloc- 
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ity field w with r = 2n/uj being the SW oscillation period. At mid- height the condition 
f(t) = —fit + t/2) implies in particular that positive and negative field extrema of an 
oscillation cycle have equal magnitudes. SWs with smaller frequency that are located in 
the end region of the SW solution branches break this symmetry (cf. Sec. IIIID|) . This is a 
prerequisite for period doubling j^. 



C. Supercritical growth and SW— >TW transition 

Here we review the supercritical spatiotemporal behavior of SW transients. We start 
from the quiescent conductive state at r < r osc , disturb it slightly by adding small random 
numbers in the range of [— 10~ 4 , 10~ 4 ], and simultaneously increase the control parameter 
r slightly above threshold as indicated in Fig. ^ Then one observes a generic transition 
scenario that occurs similarly also in analogous experiments. It consists of three evolution 
phases that we have found to be generic for convective growth out of small unspecific pertur- 
bations: (i) an exponentially growing SW of high frequency over a time interval the length 
of which depends on the size of the initial perturbations of the conductive state, (ii) an 
intermediate phase that is always very short with a spatiotemporally complicated transfor- 
mation from SW into a high-frequency TW and finally (m) a long-term TW transient to a 
low-frequency, strongly nonlinear, relaxed TW or SOC state depending on r [lfij ]. 

In generic, i.e., non specific initial perturbations that break the mirror symmetry x — ► —x 
only weakly the two critical oscillatory modes of right and left traveling waves that can grow 
above the Hopf threshold are contained with roughly equal amplitudes, i.e., \Ar\ ~ \Ai\. 
Thus, initially, i.e., as long as linear theory applies, the two critical TW modes that started 
with \Ar\ ~ \Al\ grow exponentially and independently of each other with the same growth 
rate; and their superposition causes SW-like oscillations with the large critical frequency 
as in Fig. 03 a. Hence, "almost" mirror symmetric experimental setups that do not favor a 
mrticular TW propagation direction cause an initial growth phase with SW characteristics 

421- 

But then a competition between the two TW constituents sets in when the advective 
nonlinearities have become sufficiently strong. They amplify the mirror-symmetry breaking 
differences between \Ar\ and \Al\ and cause the decay of the minority TW. Thus, the SW, 
which still has a large frequency, is transformed into a fast TW (Fig. 0b -h). This SW — ► TW 
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transformation being advection driven is spatiotemporally complicated, in particular for the 
concentration field and it implies a dramatic redistribution of concentration by advective 
"rolling in" of concentration (Fig. c -g). It takes place within less than one vertical 
thermal diffusion time and it starts roughly when the flow amplitudes of the two constituent 
TWs have grown to about the phase velocity v p : First, concentration is advected up- and 
downwards in the form of plumes by the growing SW-like flow that reverts periodically 
its direction; the vertical concentration gradient being still as large as that of the quiescent 
conductive state. The SW — > TW transformation is triggered by an advective wave breaking 
and wave toppling process of the crests and troughs of the concentration wave (Fig. c -d), 
whereas the waves of w and ST do not undergo substantial structural changes. Advective 
nonlinearities have by now become sufficiently strong to make the mirror-symmetry-breaking 
differences between the original, left and right propagating TW constituents of the SW 
clearly visible. When the concentration wave crests (troughs) with a high (low) alcohol 
content bend and topple, they are advectively "rolled in". 

The flow induced sequence of first deforming, then bending, and finally rolling in the 
plume-like wave crests and troughs is associated with and driven by a growth of the spatial 
phase shift between velocity and concentration field from zero to about A/4 during the SW 
— > TW transformation: In the SW the spatial location of the nodes of SC and w coincide 
while their oscillations are shifted in time by about a quarter of an oscillation period. The 
value of this phase difference, (pc — <p w — vr/2, does not change during the whole transition 
sequence but in the TW, it also implies a spatial shift of SC and w of A/4. 

In the early TW phase (Fig. El h) the lateral concentration difference between adjacent 
roll like regions is still almost as large as the initial vertical concentration contrast in the con- 
ductive state. Thus, the frequency and the phase velocity of this emerging TW are still very 
large, i.e., not much smaller than the critical values. But then for the parameters of Fig. El 
a long-term TW transient to a low-frequency, strongly nonlinear and strongly anharmonic 
relaxed TW state (Fig. Eli) sets in: slow diffusion degrades and homogenizes the concentra- 
tion striations, the spatial extension over which SC is constant at the two alternating high 
and low levels increases, the plateau height decreases, the width of boundary layers between 
these plateaus shrinks, the <5C-wave profile becomes more and more trapezoidal, and the 
alcohol surplus (deficiency) in the cold top (warm bottom) part of the fluid layer decreases, 
thereby reducing the overall vertical concentration difference between top and bottom. This 
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longtime degradation of concentration gradients is reflected by a dramatic decrease of the 
mixing number M and with it of the frequency u relative to the initial values — the better 
the fluid becomes mixed the smaller is to. 



D. Standing waves 



weakly nonlinear analysis (4.' 



TW and LTW convection has been studied experimentally and theoretically in detail 

141 ] . But little is known about nonlinear SW states beyond a 
j] that is restricted to the immediate vicinity of the oscillatory 
threshold. It showed that SWs are unstable there, typically bifurcating backwards. 

In order to obtain the solution branch, we stabilize the SW states by suppressing phase 
propagation (so TWs can not compete against SWs) and, if necessary, by exerting a control 
procedure. The latter operates via the field amplitudes or the heat current injected into 
the fluid and the Rayleigh number in response to the instantaneous frequency uo and its 
temporal derivative Au/At , respectively j^. In this way we trace out the SW solution 
branch all the way from close to onset with large frequency to slowly oscillating SWs that 
eventually period-double into chaos. 

Fig. El shows the bifurcation behavior of SWs and how it changes with varying Soret 
coupling strength. The solution branch for the SOC is included for comparison only for if) = 
—0.03. The heating range in which SWs exist increases when ip becomes more negative since 
the oscillatory bifurcation threshold r osc is shifted stronger to higher r than the SW saddle- 
node positions at rf w which marks the lower end of the r-interval containing SWs. All these 
SWs bifurcate subcritically out of the conductive state as unstable solutions. They become 
stable via saddle-node bifurcations when the phase-pinning condition is imposed. However, 
when this condition is lifted completely then SWs decay by developing TW transients since 
any spatial phase difference between 5C and w causes the extrema of the latter to be "pulled" 
towards the solutally shifted buoyancy extrema. Depending on r these transients either end 
in a nonlinear TW or SOC or the conductive state. 

Moving along an SW branch the maximal vertical upflow velocity w max [Fig. Ufa)] does 
not increase monotonically as in TWs and SOCs but rather has a maximum somewhat 
below the respective SOC value before it drops again. On the other hand, u and M decrease 
monotonically starting with the Hopf frequency ujh [upper end of the curves in Fig.H^b)] and 
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M — 1, respectively, at onset. This variation follows the universal scaling law uo/ujh — M, 
cf. Fig. El 

The typical bifurcation behaviour of two representative SWs of Fig. |H1 with tp=-0.03 and 
-0.25 is displayed in more detail in Fig. Full (dashed) lines refer to stable (unstable) 
branches. Note that the large-amplitude upper SOC solution branches [the one for ^=-0.03 
is shown explicitly in Fig. Ufa)] are stable down to their respective saddle-nodes when phase 
propagation is suppressed so that TW solutions do not exist. Then, the SW solutions, too, 
change their stability at saddle-node bifurcation(s), one for -^=-0.03 at r^ w pa 1.0373 and 
three for ^=-0.25 [two of them can be seen in Fig.O^b) at r^ w pa 1.1238]. With increasing r 
the flow amplitude of the stable SW slightly decreases. Then the MTS breaks. Thereafter, 
for example, the downflow (upflow) extrema occurring in the SW oscillations, say, at x=0 
(±A/2) are more intense than the upflow (downflow) extrema. This is reflected by the first 
splitting of the SW solution branches in Fig. [7| Consequently, the time averaged fields have 
now a net SOC-like structure with non zero mean downflow (upflow), say, at x=0 (±A/2). 

In Fig. [HI we show for if) = —0.03 how MTS breaking changes the SW phase dynamics 
using w, w at the mid position x = = z and M as characteristic local and global quantities, 
respectively. By definition M oscillates with twice the SW frequency as long as MTS holds 
(dash-dotted lines in Fig.[HJ). The particular MTS-broken SW orbits of Fig. |H] move closer to 
the low-amplitude fixed point on the lower SOC solution branch of Fig. Eta) with downflow 
at x — 0, so the SW spends more time in the downflow phase at x = than in the upflow 
phase. This can be seen more clearly in the temporal oscillation profile of w (Fig. 4 in [lfil ]) 
; it develops a plateau close to the value of the corresponding SOC fixed point. 

Then, at a certain r depending on if>, the first period-doubling occurs followed by further 
subsequent period-doublings that lead to chaos (Fig. [7| and Fig. |HJ). For stronger Soret 
coupling, i.e. if) = —0.25 [Fig. E£b)], we could also resolve a r-window with period-3 SW 
states and subsequent period doublings. 

Beyond the last chaotic windows, we did not observe any stable SWs: heating above this 
threshold leads to the development of transients into a stable, large-amplitude SOC on the 
■upper solution branch [e.g., with amplitude wsoc — 2.51 in Fig. Ufa), i.e., well outside the 
plot range of Fig. |H] . 

Thus, here the question arises whether - and if so - how and how long the heteroclinic 
orbits connecting the two unstable symmetry degenerate SOC fixed points on the lower 
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solution branch with small convection amplitudes (wsoc — 0.71 in Fig. |HJ) organize and 
restrain the dynamics of the periodic and chaotic SWs that switch between up- and downflow. 

IV. LOCALIZED TRAVELING WAVES 

LTWs consist of wave trains of traveling convection rolls which are surrounded by quies- 
cent fluid There exist several attemps to model LTWs by weakly nonlinear 
small-amplitude expansions around the convective onset 151, 52] and modifications thereof 

nnn 

|53L I54L 155] . But due to the strongly nonlinear characteristics of LTWs these models are 
aimed at some qualitative aspects of LTW states. 

Here we investigate stable LTWs at subcritical driving, r < r osc , where the quiescent 
conductive surrounding is stable as well. 

The rolls grow out of this environment at the tail end of the wave train, travel through 
the convective bulk of the LTW with increasing phase velocity v p (x) and wavelength A(x), 
and decay at the head of the wave train (cf. Fig. fTTT|) . These stable LTW states are uniquely 
selected. Their width I (or number of rolls) is stationary and depends in an unique way 
on the control parameters (cf. Fig. EJ). The whole convective region drifts through the 
motionless state with a drift velocity Vd that can be positive or negative but that is small 
compared to the phase velocity Q|. Therefore, the frequency of a LTW is constant in a 
frame of reference that is comoving with the LTWs drift. 

To characterize the concentration variations in a LTW we use a local mixing number 
M(x) that is defined similarly to the case of extended states as the mean variance of 5C. 
But here the lateral average is replaced by a time average at a fixed lateral position x, 

M{x) = y/(6C*)/(6(%J) . (4.1) 

So brackets denote here an average over z and t instead of z and x as in Eq. ()2.7|) . Fig.lTUlfc) 
shows the variation of M(x) over the convective region of a typical long LTW. The relation 
between M(x) and the local phase velocity v p (x) in the bulk of the LTW is the same as the 
relation between M and v p = u/k for TWs. 

The smallest stable localized pulses have for each ip a minimal width of about 5 rolls. With 
increasing r LTWs grow in width and amplitude and I diverges at a maximal driving r max . 
Beyond this value, LTWs of stationary width seem to be no longer possible. Thereafter, only 
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transient localized convective regions can occur which keep expanding between two fronts. 

As one can see in Fig. El the existence interval of LTWs at fixed if> can be divided into 
two regimes: In the regime of short pulses the frequency decreases with increasing r and 
the amplitude grows fast whereas in the regime of long LTWs the frequency increases with 
r and the amplitude grows only moderately. 

Short LTWs are dominated by a direct interaction between the growth part and the decay 
part of the convection rolls. Long LTWs are characterized by an extended bulk part where 
rolls travel with phase velocity v p (x) that increases from tail to head. This corresponds to 
an increasing local wavelength X(x) = 2ttv p (x)/uj since the frequency u is constant in the 
frame that is comoving with Vd- The speeding up of v p (x) reflects the growing concentration 
contrast between adjacent rolls along the x-direction: The minimal mixing number is located 
in the growth region of the tail where v p (x) and A(x) are minimal as well. On their way 
from tail to head the rolls do not reach a stationary balance between SC injection from the 
different boundary layers and advective mixing and diffusive homogenisation of concentration 
differences on a constant level of small M. Rather LTW rolls collapse at the head when v p 
has grown up to w max [right arrow in Fig. El (b)]. Thereafter concentration is discharged 
and sustains a barrier of (SC) ahead of the decay part which stabilizes the conductive state 
there against invasion of convection. 

Fig. El (b) and (c) compare the maximal flow velocities and the frequencies of LTWs for 
several if) with the respective TW branch of A ~ 2 for which the saddle location r T s w is the 
lowest one. Note that the LTW values for w max and uj always lie near the saddle values of 
TWs. Although there exist no extended TWs below the shown TW saddle positions there 
exist stable LTWs ahead of these global TW saddles when if) is sufficiently negative. For 
example, for if) = —0.4 almost all LTWs of finite width appear below the existence region 
for TWs as can be observed in Fig. El 

The intriguing existence of stable LTWs without coexisting TWs is ensured by a flow- 
induced lateral concentration redistribution over its convective bulk. Positive "blue" (neg- 
ative "red" ) concentration deviation from the global mean is sucked from the top (bottom) 
boundary layer into right (left) turning rolls as soon as they become nonlinear during their 
growth [left arrow in Fig. El (b)]. Convection in the bulk of a (long) LTW shows all the char- 
acteristics of a strongly nonlinear extended TW. Thus, as discussed in Sec. IIII Al positive 
(negative) concentration is transported within the closed streamline regions predominantly 
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in the upper (lower) half away from the tail towards the head. This lateral concentration 
transport is reflected in the time averaged current of 5C [cf. its streamlines in Fig. 110( e)]. At 
the same time, mean concentration, 5C ~ 0, migrates mostly to the left along open velocity 
field streamlines that meander between the closed roll regions [green regions in Fig. 110( a)]. 
In this way a region with well mixed fluid, i.e., almost vanishing 5C is created and sus- 
tained under the tail. This strongly nonlinear large scale concentration redistribution that 
maintains the concentration locally at a homogeneous mean level under the trailing front 
makes there the growth of convection possible: By reducing the stabilizing concentration 
gradients the driving buoyancy force is locally increased there to levels which suffice to cause 
convection in a well-mixed fluid. 

Fig. (d) shows that this effect causes the mean convectively generated C profile to 
extend significantly further into the conductive region than the mean convective temper- 
ature field. Thus, the buoyancy (b) [cf. Fig. [TH] (d)] is determined in the front regions 
predominantly by the concentration field. This explains (i) the decrease of buoyancy below 
conduction levels ahead of the decay part with the associated restabilization of conduction 
there and (ii) the increase of (b) out of the conductive state at the tail and its strong over- 
shoot over the bulk enabling convection growth even for heating r where no stable extended 
TWs are possible. 

V. CONCLUSION 

Three oscillatory steady states of convection in binary fluid mixtures like ethanol-water 
were investigated and compared numerically in a cross section perpendicular to the roll axes: 
Traveling waves, standing waves, and localized traveling waves. In addition the supercritical 
growth of SWs out of perturbations and the SW — > T W transformation was elucidated. 

The full bifurcation behavior of TWs and SWs has been presented and their typical con- 
centration field structures are shown in detail. Both extended oscillatory periodic patterns 
bifurcate subcritically out of the conductive state in a common Hopf bifurcation. The de- 
crease of the phase velocity or the frequency along the solution branches is connected to 
the decrease of the mean 5C contrasts in the system due to advective mixing. Thereby, the 
different symmetries and flow topologies of TWs and SWs cause a different behavior in the 
strongly nonlinear regime. However frequencies and mixing numbers are always related to 



February 8, 2008 



17 



each other via the universal scale relation uj/uj h — M. 

TWs are mirror-glide symmetric under (x, z) — > (x + A/2, — z). At high amplitudes they 
are characterized by the occurrence of closed streamlines in the frame comoving with the 
TWs phase velocity. The regions of closed streamlines are displaced alternately towards the 
upper and lower plates with their respective concentration boundary layers. The latter feed 
selectively concentration only into the nearby closed streamline domains where it is homog- 
enized. This leads to the characteristic concentration distribution with homogeneous rolls of 
alternating signs in SC. They are separated by meandering inner boundary layers of 5C ~ 
which are spatially correlated with open streamlines of the velocity field. With increasing 
amplitude the concentration contrasts between the rolls decrease, the phase velocity slows 
down, and the TW ends at large amplitudes in a mirror symmetric SOC state. 

SWs are always laterally mirror symmetric. The solution branch turns around at a saddle 
node and becomes stable there under phase pinning conditions. For strongly negative ip the 
stable part of the SW branch is directed towards lower amplitudes - in contrast to the 
TW case - and develops further saddle nodes. The complicated bifurcation behavior at low 
frequencies seems to be due to an interaction with the unstable low amplitude SOC solutions. 
There is a breaking of the mirror-timeshift symmetry which is a prerequisite for the following 
period doubling cascade into chaos. Without phase pinning the SWs develop TW transients 
which end either in a nonlinear TW or a SOC state depending on the parameters. This 
mechanism is similar to the symmetry breaking of a weakly nonlinear SW during its transient 
growth at supercritical heating towards a nonlinear TW. The transition is characterized by a 
dramatic advective reorganization of the SC field and therefore of the buoyancy field within 
a short time interval. 

The only global symmetry operation for a LTW is the mirror-timeshift [(x, z, t) — > 
(x, — z, t+r/2)] in the frame which is comoving with its small drift velocity Vd- The structual 
properties of a generic long LTW were presented as well as the global bifurcation behavior 
with r for different ip. For strongly negative separation ratios LTWs exist monostably at 
low heatings r where no extended TWs are possible. This is ensured by a flow-induced 
lateral large scale redistribution of concentration over its TW dominated convective bulk: 
Thereby, concentration is maintained under the tail of the wave train at a homogeneous 
mean level. This increases the driving buoyancy forces there locally to levels which suffice 
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FIG. 1: Bifurcation properties of maximal vertical flow velocity w max (a) and frequency u (b) versus 
r for representative relaxed nonlinear convective SOC, TW, SW, and LTW states. Full (dashed) 
lines refer to stable (unstable) states. Symbols identify the states that are discussed in more 
detail. TW and SW solutions bifurcate subcritically out of the conductive state at the oscillatory 
threshold r osc with common Hopf frequency. Vertical arrows indicate the supercritical transient 
of Fig. El from a growing SW perturbation of the quiescent conductive state into a nonlinear TW. 
The wavelength of SOC, TW, and SW is A = 2. Parameters are L = 0.01, a = 10,^ = -0.25. 
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FIG. 2: Wave number dependence of TWs at onset and at their saddle-nodes for different ijj. (a) 
Phase velocities v p (k) at the onset (straight lines) and at the saddle (dashed lines), (b) Location of 
the onset r osc (k) (straight lines) and of the saddle-nodes rj w (k) (dashed lines) for TW branches in 
the k — r-plane. The dotted line marks the onset of SOCs in pure fluids. Parameters are L = 0.01, 
o = 10. 
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FIG. 3: Mixing number M versus reduced oscillation frequencies of TWs and SWs for various 
Soret coupling ij). For SWs the time averaged mixing number is plotted. Here uh is the respective 
Hopf frequency at onset. Parameters are A = 2,L = 0.01, a = 10. 
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FIG. 4: Concentration distribution of selected extended states in vertical cross-sections through the 
rolls. TW1(2) and SOCl(2) are states marked by open (filled) circles and triangles, respectively, in 
Fig. ^ The SW refers to the lozenge in Fig. ^ The top figure in the left column shows the lateral 
wave profiles of 5C, 5T, and w at midheight, z = 0, of the fluid layer in the TW2. Parameters are 
if) = -0.25, L = 0.01, a = 10, A = 2. 
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FIG. 5: Time ordered snapshots of the concentration redistribution in the SW— >TW transforma- 
tion. The concentration distribution in a vertical cross section of the fluid layer is displayed for one 
wavelength. Full isoconcentration lines are shown for 5C > 0. Dashed lines are streamlines, i.e., 
tangents to the instantaneous velocity field. The final TW (j) propagates to the left. Parameters 
are L = 0.01, a = 10, ip = -0.25, r = 1.42, and A = 2. 
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FIG. 6: Bifurcation properties of SWs for V= -0.03, -0.05, -0.1, -0.15, -0.2, -0.25, -0.3, -0.35, -0.4 
(from left to right) : (a) Maximal vertical velocity w max . (b) Frequency u>. Unstable SWs bifurcate 
subcritically out of the quiescent conductive state [lower ends of the curves in (a); upper ends in 
(b)] and undergo stability changes via saddle-node bifurcations. The SOC solution branch is shown 
for the sake of clarity only for ifj= -0.03. SOC curves for the other tp are shifted slightly to the 
right. Parameters are L = 0.01, a = 10, and A = 2. 
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FIG. 7: Details of the SW bifurcation behavior for ip = —0.03 (a) and ip = —0.25 (b). The 
magnitude |w ea: fr| of the extrema in the vertical flow is shown versus the reduced Rayleigh number 
r. The SWs bifurcate subcritically out of the quiescent conductive state at r osc . When phase- 
pinning conditions are applied they become stable (solid lines) at saddle-nodes positions [lying 
outside the plot range of (b)]. When the MTS is broken, the solid SW line splits into two since 
the magnitudes of the vertical flow extrema occurring during one oscillation cycle become different 
[see, e.g., Fig. |S] (b) where the downflow at x = = z is more intense than the upflow]. This 
MTS-broken SW starts at the arrows to undergo a period-doubling scenario (cf. insets) leading to 
chaos. Parameters are L = 0.01, a = 10, and A = 2. 
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FIG. 8: MTS-breaking and period doubling in the phase space dynamics of SWs. Shown are w 
and the mixing number M versus the vertical velocity w at x = = z for -0=-O.O3. In (a)-(g) 
the full (dash-dotted) lines refer to asymmetric (symmetric) SWs for which the MTS is (not yet) 
broken. Period doubling is displayed in (b)-(d) and in (f)-(g). A chaotic trajectory is shown in 
(h). Upwards and downwards pointing triangles indicate symmetry degenerate unstable SOC fixed 
points with upflow and downflow, respectively, at x = that lie on the low-amplitude SOC solution 
branch [dashed line in Fig. [Tfa)]. The flow velocity on the stable large-amplitude SOC solution 
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FIG. 9: Bifurcation properties of LTWs (symbols) and TWs (lines) for different separation ratios 
ip: (a) Full width I of LTWs at half maximum of the envelope of the vertical velocity field w. (b) 
Maximal vertical flow velocity w max . (c) Frequency (for LTWs in the frame comoving with the 
drift velocity Vd of the respective LTW). Filled circles identify the LTW whose structure is shown 
in Fig. I1UI Lines in (b, c) denote TWs with saddle-node wave number kj w ~ n. Unstable TWs 
(dashed lines; determined with a control method) bifurcate subcritically with large Hopf frequency 
uji at r osc out of the conductive state and become stable (solid lines) at the saddle-node rj w when 
lateral periodicity is imposed with A = 2ir/k. Parameters are L = 0.01, a = 10. 
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FIG. 10: Broad LTW of l=YI A identified by filled circles in Fig. |HJ (a) Snapshot of concentration 
deviation 5C from the global mean value in a vertical cross section of the layer, (b) Snapshots of 
lateral wave profiles at midheight, z=0, of SC, vertical velocity w, and its envelope. At the arrows 
Wmax = Vp. (b) Mixing number M(x), [Eq. (|4.1|l ] and phase velocity v p of nodes of w(z = 0) in the 
frame comoving with the small drift velocity of the LTW. The variation of \(x) = 2ir v p (x)/uj is 
the same since the LTW frequency u is a global constant, (d) Time averaged deviations from the 
conductive state at z=-0.25 for concentration, temperature, and their sum, b = T—T con d+C—C con( i, 
measuring the convective contribution to the buoyancy, (e) Streamlines of time averaged velocity 
field (u) and concentration current (J) = (uSC — LV(5C — tpST)). (u) results from (b) and affects 
(J) via the contribution (u)(SC). In the upper half of the layer positive 8C (alcohol surplus) 
is transported to the right. In the lower half of the layer negative SC (water surplus) is also 
transported to the right - positive 5C is transported there to the left as indicated by the arrow. 



